Theory of Current-Induced Magnetization Precession 



M.D. Stiles 

National Institute of Standards and Technology, Gaithersburg, MD 20899-8418 



o 
o 

(N 



C/3 



Jiang Xiao and A. Zangwill 
School of Physics, Georgia Institute of Technology, Atlanta, GA 30332-0430 
(Dated: February 2, 2008) 

We solve appropriate drift-diffusion and Landau-Lifshitz-Gilbert equations to demonstrate that 
unpolarized current flow from a non-magnet into a ferromagnet can produce a precession-type 
instability of the magnetization. The fundamental origin of the instability is the difference in 
conductivity between majority spins and minority spins in the ferromagnet. This leads to spin 
accumulation and spin currents that carry angular momentum across the interface. The component 
of this angular momentum perpendicular to the magnetization drives precessional motion that is 
opposed by Gilbert damping. Neglecting magnetic anisotropy and magnetostatics, our approximate 
analytic and exact numerical solutions using realistic values for the material parameters show (for 
both semi-infinite and thin film geometries) that a linear instability occurs when both the current 
density and the excitation wave vector parallel to the interface are neither too small nor too large. 
For many aspects of the problem, the variation of the magnetization in the direction of the current 
fiows makes an important contribution. 
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I. INTRODUCTION 

Recently, the phenomenon of giant magnetore- 
sistance has been used to unambiguously detect 
current-induced magnetization reversal in thin magnetic 
heterostructuresiSi^ These observations confirm theo- 
retical predictions^ that angular momentum can be 
transferred from current-carrying electrons to the mag- 
netization of a ferromagnetic film. Earlier reports^ of 
the observation of this so-called spin-transfer effect at- 
tributed a magnetic origin (magnetization reversal, free 
precession, or spin wave excitation) to anomalies seen in 
transport measurements. All these experiments used a 
ferromagnet/non-magnet/ferromagnet spin valve geome- 
try because a "polarizer" ferromagnet was deemed nec- 
essary to ensure that the electrons incident on the "an- 
alyzer" ferromagnet from the non-magnet would carry a 
spin-polarized current. Now, two new experimental facts 
have appeared that bear directly on these issues. First, 
current injection into a single ferromagnetic layer has 
been observed to produce resistance anomalies very sim- 
ilar to those seen in the spin valvesi Second, microwave 
emission has been detected in at least some spin valve 
samples when spin-transfer resistance anomalies occur^ 

Based on the foregoing, it is natural to inquire whether 
magnetic precession occurs when an unpolarized cur- 
rent flows (perpendicular to the interface) from a non- 
magnetic metal into a ferromagnetic metal. We investi- 
gate this question here by performing a linear stability 
analysis of a presumptive current-induced precessional 
state. Our approach to this problem of coupled trans- 
port and magnetization dynamics derives from previously 
published work by two of us. One earlier papes^ dis- 
cussed the quantum mechanical origin of spin-transfer 
torque at a ferromagnetic/non-magnctic interface. An- 
other papesifi solved a matrix Boltzmann equation to 
compute spin currents, spin accumulation, magnetoresis- 



tance, and spin-transfer torques for multilayers. In the 
present work, we study both single and multiple interface 
geometries in the presence of an external magnetic field 
(normal to the interface) strong enough to saturate the 
magnetization (Figure^ . Magnetic anisotropy and mag- 
netostatics are neglected for simplicity. Our approximate 
analytic and exact numerical solutions reveal that an ini- 
tial instability toward precession occurs when both the 
current density and the excitation wave vector parallel 
to the interface are neither too large nor too small. Oth- 
erwise, Gilbert damping suppresses the instability. An 
important feature of our model is that the magnetiza- 
tion is permitted to vary in the direction normal to the 
interface. Among other things, this means that preces- 
sion can occur in the film geometry even when the two 
interfaces are identical. Other systematic features of the 
precessional state will be described in detail below. 

Our results complement and extend previous theo- 
retical studies of current-induced precession in mag- 
netic films. In their pioneering work, Slonczewski^ 
and Berger^- used analytic model calculations to pre- 
dict the critical current for the onset of precession in 
spin valve structures with uniform but non-coUinear mag- 
netizations in the two ferromagnetic layers. Recent 
work shows that the seemingly different answers they 
obtained are not inconsistent lii Two global studies of 
precession and switching in single-domain spin valvesi^ 
and single-domain particlesi^ supplement the Landau- 
Lifshitz-Gilbert equation with the spin-transfer torque 
derived in Ref. 0. Bazaliy et al.— relaxed the assump- 
tion of uniform magnetization and studied the effect of a 
legislated spin current flowing from a non-magnet into a 
semi- infinite ferromagnet. They found that spontaneous 
excitation of spin waves occurs when the current den- 
sity is sufficiently large. Finally, we draw attention to 
a very recent application^^ of magnetoelectronic circuit 
theorjii^ to study current-induced precession in an ultra- 



Jo 



- H 



(a) 



Non-magnet ^ Ferromag^^^J 


(b) 


t 

< ► 












t 

M ► 











FIG. 1: The external magnetic field and particle current both 
point in the positive z-direction for (a) a single-interface ge- 
ometry; (b) a thin film geometry; and (c) an electron reser- 
voir geometry. The interfacial areas are infinite in all cases. 
In (b), the non- magnetic leads are semi-infinite in length. In 
(c), the non- magnetic leads are attached to electron reservoirs 
at equilibrium. 



thin ferromagnetic film with uniform magnetization. For 
unpolarized current densities above a critical value, the 
authors predict that unstable spin wave modes exist only 
for non-zero wave vectors parallel to the interface. Our 
results confirm this prediction. 

The calculations presented in this paper combine phe- 
nomenological constitutive equations for the electric cur- 
rent and the spin current (drift-diffusion equations) with 
a Landau-Lifshitz- Gilbert equation generalized to include 
spin-transfer torque. The boundary conditions at each 
non-magnet/ferromagnet interface are chosen so that rel- 
evant observables reproduce the results of a Boltzmann 
equation treatment of the same problem. Using experi- 
mental data to fix the parameters, we choose a value for 
the electric current density, and solve for the interface 
charge density (chemical potential discontinuity), spin 
accumulation, and spin current assuming a precessing 
magnetization with an amplitude and phase that depends 
on the running variable z in Figure^ The presumed pre- 
cession is a legitimate mode of the system if the computed 
spin-transfer torque balances the ferromagnetic exchange 
torque at each interfacei^^ The mode is unstable (spon- 
taneous excitation in the presence of a current) if the 
imaginary part of the mode frequency is negative. The 
mode is stable and uninteresting if Gilbert damping over- 
whelms the effect of spin transfer. 

The plan of this paper is as follows. In Section II, 
we set down the basic definitions and equations of the 
theory. Section III gives a qualitative picture of the origin 
of the precessional instability. The boundary conditions 
for the transport and magnetization equations are given 
in Section IV. Sections V and VI present the details of the 
calculation for a single interface geometry and a thin film 



geometry, respectively. Section VII discusses our results. 
Section VIII examines the effect of non-magnetic leads 
using a reservoir geometryii^ Section IX discusses some 
limitations of our model and suggestions for future work. 
Section X is a summary. Appendix A derives the effective 
exchange field for a bounded ferromagnet. The transport 
boundary equations are derived in Appendix B. 



II. DEFINITIONS & BASIC EQUATIONS 

In an earlier paper we defined the fundamental vari- 
ables of spin transport quantum mechanically and de- 
scribed the various relations that exist among them. 
Here, we proceed phcnomcnologically and decompose the 
total number density of conduction electrons into major- 
ity and minority contributions n = n-\ + ni so that the 
current-induced deviations of these quantities from their 
equilibrium values are related by 



5n = 5n^ + 6n\^. 



(1) 



In the drift-diffusion approximation, the particle current 
in the non-magnet is 
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where E is the electric field, a is the conductivity, and 
D is the electron diffusion constant. We speak of spin 
accumulation in the non-magnet when the spin density 



(3) 



is different from zero. Notice that we use the particle 
current and the spin density as variables rather than the 
charge current and the magnetization. 
A so-called spin current density 
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flows in the non-magnet when the spin accumulation is 
not uniform in space. In Q), the first index of 
labels the a**^ Cartesian component of the electron spin 
vector while the second index labels the /J**^ Cartesian 
component of the flow direction. 

In the ferromagnet^ we make explicit use of the 
fact that there is no appreciable spin accumulation 
or spin current in the direction perpendicular to the 
magnetizationAii Accordingly, 5m ~ {h/2){5n-\ — Sn^) 
is the spin accumulation so 



M = {nis + dm)u 



(5) 



is the total spin density. Here, u is a unit vector, mg — 
hMs/gfiBi a-iid is the saturation magnetization. We 
adopt a two-current model^^ so j^'^ = j| + j| where 

j| (cr|/e)E - Z^iV^ni 



JX = iai/e)E-D_iVSn 



(6) 
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using an obvious notation for the conductivity and dif- 
fusion constants of majority and minority spin electrons. 
Combining all the above, the total particle current in the 
ferromagnet is 

•i'M^_\ \ iV(5n '— — ^V5m. (7) 

The corresponding spin current density in the ferro- 
magnet is 

Qlf^h^n^ih-h)p- (8) 

Consistent with our neglect of transverse spin accumula- 
tion in 0, this formula presumes that the spin part of 
the spin current in a ferromagnet is always parallel to the 
spin density. Then, using @, we get 



where fsf is the spin-flip scattering time in the ferromag- 
net. The effective field Hgff accounts for the external 
field, exchange, and damping: 
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(13) 



Here, lex — ■\/2Acx//ioA^| is the exchange length in a 
ferromagnet with spin stiffness Acx, wo — T^fs, and a 
is a phenomenological damping parameter. We ignore 
magnetostatics, but for fields and magnetization perpen- 
dicular to the interface, its main effect is to shift the 
applied magnetic field by an amount Mg. 



III. ORIGIN OF THE INSTABILITY 
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Textbooks show that the charge density inside a 
current-carrying Ohmic conductor is zero (except within 
a screening distance of an interface or free surface). 
Therefore, we are justified to put (5n = in |5J, 0, and 
0. We will also fix the current to be j = j^z everywhere. 
Using this fact in Q, we can eliminate E from © and 
simplify the components of the ferromagnetic spin cur- 
rent density that flow in directions longitudinal (z) and 
transverse (±) to the interface normal. The result is 
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Two characteristic features of the spin current Qa/j 
underlie the qualitative physics of current-induced pre- 
cession. First, the structure of the spin-transfer torque 
in Hll|) and H12|l implies that the magnetization at the 
z = interface in Figure ^ feels a torque per unit area. 
It is 
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lim I dzV Q = (Q^M - Q™) 
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(14) 



where ±' refers to the directions transverse to the fer- 
romagnetic magnetization (z' is the direction along the 
ferromagnetic magnetization). The last line in (|14|l is 
true because 
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(10) 



where 7P = {a^ — (c]- + cr^ ) is the degree of cur- 
rent polarization desired by the ferromagnetic bulk, — 
(D|(T4 + Dia^)/{a^ is an effective longitudinal spin 

diffusion constant, and D± = |(I?| + -Dj) is an effective 
transverse spin diffusion constant. 

The equation of motion for the spin density in 
the non-magnet is a continuity equation for the spin 
accumulation^ Taking account of distributed spin- 
transfer torque, spin flip scattering and torque from the 
external magnetic field, we find 
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-7[mxHe.tL + — = 0, (11) 



where Tgt is the spin-flip scattering time in the non- 
magnet and 7H = {gfiB/h)'B. The equation of motion 
for the spin density in the ferromagnet is the Landau- 
Lifshitz- Gilbert equation supplemented by a term to ac- 
count for spin-transfer torque as in (|11() . Specifically, 
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and Q™ ~ at the interface between a non-magnet 
and a ferromagnet due to spin-filtering and other 

The second bit of information we need is a linear ki- 
netic equation that relates the transverse spin current 
at the interface, Q^5^, to the interfacial densities of our 
problem (charge accumulation, spin accumulation in the 
non-magnet, and ferromagnetic spin accumulation). Of 
these, for a static magnetization only the spin accumula- 
tion m in © has a transverse component. Therefore, it 
must be the case that 



QNM _ L ™ 
— omj_/ 



(16) 



where 5 is a constant. More detailed considerations (see 
Appendix B) show that 6 > 0. 

For simplicity, we will study (in this section) a time- 
independent situation where electrons flow from a non- 
magnet into a ferromagnet with a small amplitude, frozen 
spin wave at the interface. That is, in (jSJ, we let u = 
z -|- exp(iK • R) where R and K are two-dimensional 
vectors in the plane of the interface. This means that 
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FIG. 2: (a) Cartoon of steady-state spin configuration when 
K — 0; (b) spin accumulation and spin current when K = Q; 
(c) sohd curves: mz{x = 0, z) and mx{x = 0, z) when K — 0; 
dashed curve: mx{x = 0,2) when K ^ 0; (d) sohd curves: 
Q^f{x = 0,z) and Q^fix = 0,2) when K = 0; dashed 
curve: Q^f (a; = 0,z) when AT / 0; (e) same as (b) when 
K 0; (f) same as (a) when K ^ 0. See text for details. 



conduction electrons incident on the ferromagnet from 
the non-magnet see a static magnetization 



M = M{ReKe*^-^],ImKe*^-^],l} 



(17) 



In the non-magnet itself, l@J and Ijlll) with Hcxt — 
combine to give 



Tsf 



(18) 



for the steady-state spin accumulation in the non- 
magnet. This is solved by 



m 



= {Rc[mi(z)e'K-^],Im[mi(z)e'K.R]^^^(^)|^ (^g) 



if 



mz{z) = mzexp(z/Zsf) 
m^(z) = m^cxp(Kz) 



(20) 



with = DTsi and Z^jK = 1 



We begin our analysis with the right side of the top 
panel in Figure [3 This shows a spin pattern in the fer- 
romagnet that is spatially uniform {K = 0). In panel 

(b) , just below, the arrow labeled m reiterates the infor- 
mation conveyed by the diffusing spins on the left side 
of panel (a). Namely, the spin accumulation in the non- 
magnet is everywhere anti-parallel to M. To understand 
this, we note first that Q™ > on account of the first 
term in p(J|) .^° Therefore, to satisfy p5|) . there must be 
a collection of spins in the non-magnet (polarized paral- 
lel to M) flowing from left to right toward the interface. 
We interpret this, equivalently and consistently with the 
diffusive solutions H2U|) graphed as solid curves in panel 

(c) , as a flow of anti-parallel spins away from the inter- 
face into the non-magnet. Qualitatively, we imagine that 
spins diffuse into the non-magnet after being "imprinted" 
in the anti-parallel direction by the ferromagnet. Ran- 
domly walking spins (including those that happen to dif- 
fuse back to the interface) only encounter spins parallel 
to the themselves because all imprinting was done by fer- 
romagnetic spins that are parallel. From (|14|l and 

Nst — AQ^^ = because m_L' = 0. There is no driving 
force for precession when the ferromagnetic magnetiza- 
tion is uniform. On the other hand, the spin currents 
and are both non-zero [solid curves in panel 

(d) ] because M is not normal to the interface. 

The situation changes markedly when the magnetiza- 
tion M varies along the interface {K ^ 0) [right side of 
Figure[2If)]. before, spins diffuse into the non-magnet 
after being imprinted in the direction anti-parallel to the 
local ferromagnetic magnetization. However, because the 
imprinting ferromagnetic spins now differ in their trans- 
verse directions, spins diffusing in the non-magnet typi- 
cally encounter spins whose transverse component differs 
from their own. The net effect is a reduction in the trans- 
verse spin accumulation [dashed curve in panel(c)]. The 
effect is greater farther away from the interface because 
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FIG. 3: (a) Interfacial spin current and spin for current flow 
from the non-magnet into the ferromagnet. This diagram is 
the same as Figure |2Ie) ; (b) Same as (a) except that the 
current flows from the ferromagnet into the non-magnet. 



the transverse diffusion is superimposed on longitudinal 
diffusion away from the interface. From Q), there is 
a corresponding increase in the transverse spin current 
near the interface [dashed curve in panel (d)]. Compared 
to panel (b), the spin accumulation vector m in panel 
(e) bends closer the negative z-axis and thus acquires 
a positive component perpendicular to M. Therefore, 
from H14() and (|16|) . the vectors and Ngt appear as 

sketched in panel (e). The indicated torque tends to ro- 
tate M away from the interface normal. That is, Ngt 
tends to increase the amplitude u of the incipient spin 
wave assumed at the beginning of the discussion. This 
is the fundamental origin of the precessional instability. 
Thus, in contrast to the common spin- valve geometry, we 
have here a single ferromagnet that both generates a spin 
current and is acted on by that spin current. 

Two additional qualitative points emerge form the dis- 
cussion above. First, if the current jo reverses direc- 
tion, so electrons flow from the ferromagnet into the non- 
magnet, H10|l tells us that the sign of Qzz reversest22. The 
direction of the spin accumulation reverses also so Q 
remains correct. The decrease in transverse spin accu- 
mulation which rotates m toward the z-axis now implies 
that the direction of Ngt reverses as well (Figurel^)). This 
torque tends to reduce the amplitude of the presumed 
spin wave. We conclude that electron flow from a fer- 
romagnet into a non-magnet tends to suppress a preces- 
sional instability. This type of non-reciprocal behavior 
with respect to the direction of current flow is a charac- 
teristic feature of spin-transfer phenomena.'* 

Finally, we remark that the presence of a second "po- 
larizer" ferromagnetic film upstream in Figure H tends 
to enhance the precessional instability. Even the K = 
mode can become unstable as long as the magnetization 
of the polarizer is not collinear with the magnetization 
of the "analyzer" . This is so because, in that case, the 
electrons incident on the z = plane carry a spin current 
with a non-zero component in a direction perpendicular 
to M. As a practical matter, this means that the criti- 



cal current for the onset of precession is lower when two 
ferromagnets are present than for the case of a single 
ferromagnet studied by us. 

IV. BOUNDARY CONDITIONS 

A quantitative solution of our problem requires bound- 
ary conditions for the transport and magnetization equa- 
tions set down in Section II. One of these arises because 
the exchange interaction is effective right up to the sur- 
face of a bounded ferromagnet. For the thin film ge- 
ometries of Figure ^ we show in Appendix A that an 
exchange torque density riox = — 7M x exists where 

Ul^='^[5{z)-5{z-t)]dM. (21) 

This shows that a localized torque acts on each inter- 
face in Figure ^b). For the single interface geome- 
try, the 6(z — t) term is absent. Spin wave calculations 
for bounded ferromagnets^* eliminate the surface torque 
produced by H21() by imposing the boundary condition 
dzM. = at each interface. For our problem, the spin- 
transfer torque 11411 also acts at each interface. This dic- 
tates that we determine the spin wave modes by impos- 
ing a zero net torque requirement at each interfaced The 
specific boundary condition is 

Ne, + N,t = 0. (22) 

where Ncx is the exchange torque density iicx integrated 
over an infinitesimally thin pillbox as in H14() . 

A distribution of charge/spin current where the spins 
precess at a complex frequency uj = -I- iuji is a mode 
of the system if (lll|l is satisfied in non-magnetic regions, 
(|12|l is satisfied in ferromagnetic regions, and 1)221) is sat- 
isfied at every interface. When uii > 0, the mode is sta- 
ble and uninteresting because its amplitude never grows 
large enough to be observed. True current-induced pre- 
cession corresponds to modes with cui < 0. 

Of course, the currents must also be consistent with 
the spin-dependent scattering that occurs at the interface 
between a magnet and a non-magnet. Previously, two of 
the present authorsiS addressed this issue by imposing 
scattering boundary conditions on a 2 x 2 matrix Boltz- 
mann equation and computing all observables from the 
semi-classical distribution functions (k, r) and (k, r) . 
In this paper, we apply boundary conditions to the ob- 
servables directly in a manner that is consistent with this 
Boltzmann equation treatment. 

We show in Appendix B that the boundary conditions 
obtained in this way take the form of a set of linear kinetic 
equations which relate the interfacial currents to the in- 
terfacial densities. Here, we quote a form that is valid for 
small amplitude precessional motion around Hcxt — Hz. 
That is, we assume (here and in all subsequent work) 
that u = z -I- u± with u± ^ 1: 

jo = Sno[wo{X^ + Xi)] + 2h-'mz[wo{X^ - X^)] 
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SO from (jSj), 

M = M{z){Rc[u±{z)<S>],lm[u±{z)<S>],l} (28) 

where 

M{z) ^ms + 5m{z). (29) 



Q^M = iMnoK(^T-^i)]+"^.K(^T+^i)] 



The constant 



Wo 
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2 9^ 



(23) 



(24) 



where Aps is the area of the Fermi surface projected onto 
the interface. The constants and are defined sim- 
ilarly with the densities of states per spin dnij /dfi^ and 
dni/dfii- Finally, 



(25) 



where R| = I^tI^ ^^"^ ~ l^iP ^''"s the squares of 
the majority and minority electron reflection amplitudes 
averaged over the Fermi surface. The appearance of the 
interfacial charge accumulation Suq on the right side of 
(|23|1 is not inconsistent with our earlier statement that 
j^NM _ ^j^FM = Q ^j^Q bulk. There is a resistance 
associated with the interface whenever Ro- ^ 0. The 
corresponding voltage drop may be thought of as arising 
from a dipole layer of charge difference Shq localized at 
the interface. 

We note in closing that if uj^ were zero (M || z), the 
second and third boundary conditions in (|23|) reduce to 
(|16|l . Appendix B demonstrates that the terms propor- 
tional to in 1)23(1 account for the fact that M is slightly 
tilted away from the interface normal. 



V. SINGLE INTERFACE GEOMETRY 

We showed in Section III that spin-transfer torque 
can induce a precessional instability. To discover if it 
actually does so requires a full solution of the coupled 
transport/magnetization dynamics equations with their 
boundary conditions. In particular, we need convenient 
formulae for the interfacial torques that enter (|22|) . For 
the single interface geometry of Figure ^a), a suitable 
precessional ansatz for the magnetization in the non- 
magnet is 



m — {Re[m±{z)^],lm[m±{z)^],niz{z)} 



(26) 



where $ = gi(K R+(^t)_ rj,]^^ corresponding ansatz in the 
ferromagnet is 



A. Interfacial Spin- Transfer Torque 

To find the spin-transfer torque, we substitute 1)26(1 and 
into the longitudinal component of the equations of 
motion and The solutions are [cf. 



and 



with l?f 



mz{z) = m^c-xpiz/Ui) 



Sm(z) = Smexp(—z/lsf) 



(30) 



(31) 



— DzTgf. These are used in the longitudinal 
component of I^J and ((10(1 to get the interfacial relations 



■' ZZ 



-iD/ki)m, 



and 



' ZZ 



(32) 



(33) 



The previous two equations, together with the first and 
last (longitudinal) boundary conditions in l(23() . deter- 
mine the four interfacial quantities Suq, rUz, Sm, and Qzz 
in terms of jo, l"^, D/l^i, Dz/lsi, X-^, X^, wq, wt, and 
wi . In particular, all four are proportional to jo, and thus 
change sign when the current reverses as argued earlier. 

Returning to 126() . we choose Hoxt = Hz, so the trans- 
verse component of ((11() is solved by 

mj_{z) = m±exp{Kz), (34) 

where k satisfies 

rather than the formulae quoted after ((20(1 . Using 1(34(1 . 
the transverse component of yields 

Q^f ^ -Dnm^ (36) 

while the transverse boundary conditions in l(23() are 

Qff ~ u^Q^f = woim^ - u^mz). (37) 

From these, we can solve for the four interfacial quan- 
tities Q,^^ and m^. The spin-transfer torque follows 
straightforwardly from ((14(1 when we follow Appendix B 
and express Q^5^ in unprimed coordinates (z-axis paral- 
lel to the interface normal) rather than primed coordi- 
nates (z'-axis parallel to M): 

— - H±z - ^-LQzz - ^ ^ A ^ ' 



u = z -I- u^(z)$. 



(27) The constant A = wo^sf/-D. 
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B. Interfacial Exchange Torque 

The torque produced by the surface exchange field l|21|) 
over an area A at z = is computed by integrating the 
torque density — 7M x Hf^ over an infinitesimaUy thin 
piUbox as in H14(l . The result is 



Nex 

A 

because 



^^^MxdM=~i'^M\z)n\ (39) 



Z X Ux = iU± 



(40) 



corresponds to a 90° rotation of the vector uj^ . Since the 
function M{z) in (|29|l is known at this point, the only 
unknown in (|39|l is u'^ = dzU±{z)\z=o- 

We determine u'j_ numerically by integrating the 
transverse component of the Landau-Lifshitz-Gilbert 
equation H12|l inward from z = 00 to z = 0. The equation 
of interest is 







{iuj 



1 



+ ^ 



M{z 
oJoM{z)l 



ijH} uj 



{lnjPjo-DzM'{z)}u'j 



The first term in brackets comes from the Ma, damping, 
and external field terms. The second term in brackets 
comes from the spin-transfer torque in the bulk of the 
ferromagnet due to the spatial variation of the magneti- 
zation. The third term in brackets is the bulk exchange 
torque. On account of l|31() . the spin flip term u±6m/fsf 
in (|12|l cancels a piece of the spin-transfer torque. To 
derive (|^ . we have used repeatedly and exploited 
the identity 



\7^{Mu) ^ uV^M + 2VM ■ Vu -I- AfV^u. 



(42) 



Very far from the interface, M{z) 
duces to 



nis and (|^ re- 



— {iLU + aoj — i"fH — IujqK ^ex)"-L 



+ \n-i^{jQ/m,)xi\ + iLOoll^n'\. (43) 

This equation is solved by Uj^(z) = Uj^exp(— pz) where 
(Rep > 0), 



pi 



ox - -i^ ± ■^^J-n^^ -^ojo^, (44) 



2ujQ 2wo 

and A = cj — iuja — — ujoK^l^^, where 



2L 



(45) 



This solution is used to start the inward integration of 
(I41II . We note in passing that if a = = and q = ip, 
(03 simplifies to 



uj{q) = jH + uJol^^q^ 



(46) 



This is (essentially) the dispersion relation derived by 
Bazaliy and co- workers!'* for bulk spin waves in the pres- 
ence of a specified spin current. An instability occurs 
(uj goes negative) when the current is large enough that 
> Aujo-fH. If we take iJ,oH = 1 T, this gives a charge 
current density Jq = \e\jo > 10 x 10^ A/cm^, which is an 
order of magnitude larger than the interfacial instabilities 
we discuss below. 

It is now straightforward to search for precession-type 
instabilities. For given values of current jo, external mag- 
netic field iJ, and transverse wave vector K, we search 
through complex frequency {uj = ojr -|- iui) space looking 
for modes. That is, situations where H22|l is satisfied us- 
ing (|38|l for the interfacial spin-transfer torque and H39() 
for the interfacial exchange torque. If uji < 0, the mode 
is linearly unstable against precession. A suitable initial 
guess is the zero-current surface spin wave frequency 



1 — ia 



(47) 



VI. THIN FILM GEOMETRY 



As in the previous section, the search for current- 
induced precessional instabilities in the thin film geome- 
try of Figure ^b) begins by solving the longitudinal part 
of the problem. The spin accumulation decays exponen- 
tially into the non-magnetic leads so at the left (z = 0) 
and right (z — t) interfaces, we have the analog of H32() 
for the interfacial spin currents: 



£1 

If, 



(48) 



The superscript on the diffusion constants and on the 
spin fiip lengths allows us to assign different material 
properties to the leads and thus to the interfacesi^^ 

The z-dependences of the spin accumulation and the 
spin current in the ferromagnetic film are linear combina- 
tions of exp(±z/rsf). Straightforward algebra yields the 
interfacial relations 



D 



[Sni^ coth(t/Zsf) - (5m^csch(i/4t)] 



[Sm^ coth{t/lsi) + 67n^csch{t/lsf)] . 

(49) 

The four equations in H48|l and (|49|l . together with the 
two longitudinal boundary conditions in l|23|) at each in- 
terface, determine the eight quantities Suq, Sn^, , ttl^, 
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Siri^, 6m^, Q^zJ ^^'^ Q^z terms of the physical param- 
eters of the modeL FinaUy, because it is used in the 
integration of 141(1 . we record that 

M{z) = TOs + (5m'"csch(t/rsf) sinh[(i — 2)/rsf] 

+ (5TO^csch(i/Zsf ) sinh(z/4f). (50) 

Up to an overall sign, the spin transfer and exchange 
torques are still given by (|^ and l(^ . Therefore, at the 
right interface, we choose arbitrarily and compute u'j_ 
so that H22|l is satisfied at that interface. This provides 
the initial condition for the integration of 1(4111 toward z = 
0. Then, the search through complex frequency space for 
solutions that satisfy ((22(1 at the left interface proceeds 
as described earlier. As initial guesses, we use the zero- 
current frequencies of the spin wave modes of the film: 

— : • (,0-1- j 

1 — la 

VII. RESULTS & DISCUSSION 

The numerical results presented in this section are 
based on converged solutions of the transport and mag- 
netization equations using material parameters derived 
directly from experiment or from first-principle electronic 
structure calculations (Table QJ. 

A. Single Interface Geometry 

The top panel of Figure ^ illustrates the regions of 
linear instability as a function of applied charge current 
Jo and wave vector K parallel to the interface for the 
single interface geometry (top panel of Figure . An 
initial state of current induced precession occurs inside 
each contour labeled by a choice of external magnetic 
field H. We emphasize that every contour forms a closed 
loop in the Jq — K plane. That is, for a fixed value of Jq, 
the instability region is bounded by a lower and an upper 
critical wave vector and, for a fixed value of K, the insta- 
bility region is bounded by a lower and an upper critical 
current density. Note also that the instability occurs only 
for rather small values of external magnetic field. We will 
see in the next subsection that thin films support preces- 
sion for much larger values of H. In addition, it appears 
that the upper critical field is very sensitive to the param- 
eters that describe the system. A systematic exploration 
of parameter space is prohibitive, but we note that the 
choices Isi = 1000 nm and Igf = 15 nm give regions of 
instability for iJ > 2 T. 

The bottom panel of Figure 01 is the same as the top 
panel except that we have replaced the full transverse 
Landau-Lifshitz- Gilbert equation ((41(1 by 

— {iuj + auj — ijH} u± 



Parameter 


Value 


Ref. 


p(Cu) 


6 




22 


PT 


40.5 


nfim 


22 


Pi 


109.5 


nrim 


22 


kc{Cn) 


450 


nm 


23 


Lf(Co) 


59 


nm 


24 


AR^ 


0.117 


f^m^ 


22 


ARi 


0.903 


f^m^ 


22 


^FS 


580 


nm~'^ 


25 


dn/dfj. 


6.23 


eV-^nm-3 


26 




8.17 


eV"^nm~^ 


26 


dni/dfj,i 


31.6 


eV-^nm"^ 


26 




2 X 10~" 


J/m 


27 


M, 


14.46 X 10^ 


A/m 


28 


7(Co) 


1.9x10" 




29 


7(Cu) 


1.75x10" 




25 


a 


0.01 




30 


D 


4.1 X 10^^ 


nm'^s"^ 






1.7 X 10^^ 


2 —1 

nm s 




rsf(Cu) 


2.4 X 10"" 


s 




rsf(Co) 


2.0 X 10"^^ 


s 




7^ 


0.46 








3.465 xlO" 


s-i 






3.9 


nm 




2ms/7i 


143.6 


nm~^ 




Wo 


11.4 X 10" 


nms~^ 




u;| 


8.7 X 10" 


nms~^ 




Wl 


2.25 X 10" 


nms"^ 






1.90 






Xi 


0.24 







TABLE I: Material parameters for Cu/Co. The values above 
the horizontal line were taken from experiment or first- 
principles calculation. The values below the line were cal- 
culated using equations given in the text, a = 1/p, and the 
Einstein relation e^a = Ddn/dfi. 



+ —{lnjPjo~DzM'{z)}u'j_ 

+ zu;o4{uT -^^U-l} • (52) 

This approximation neglects the ferromagnetic spin ac- 
cumulation 6'm{z) in the exchange interaction and re- 
tains only its gradient variation in the bulk spin-transfer 
torque. Evidently, these omissions make no difference to 
the qualitative features of the diagram. 

According to Figure 0] no precession occurs if either 
the current density, the magnetic field, or the parallel 
wave vector gets too large. We expect the instability 
to disappear when the Gilbert damping overwhelms the 
effect of the spin transfer torque. Therefore, since the 
damping parameter a enters the theory only in the com- 
bination au [see ((52ll ]. we anticipate that the mode fre- 
quency is an increasing function of j, H, and K. 
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FIG. 4: Regions of instability for the single interface geom- 
etry of Figure ^a). Each contour gives the region in Jq-K 
space within which instabilities occur for an series of applied 
fields (labeled on each contour). The top panel gives the re- 
sults for the full calculation in the ferromagnet. The bottom 
panel gives the results if the effect of spin accumulation on 
the exchange interaction is ignored. 



To check this, we approximate the spin current in the 
second line of H52|l by its value at the interface to get 



= {iuj + acu — ijH} 



TTLs 



(53) 



This linear equation, like H43|) . is solved by uj^(z) 



u±exp{—pz). Therefore, 

with p given by (|44() but with H45() replaced by 
Wo/2 - D,SmiO)/kf h-/'Pjo 



(54) 



n = 



2L 



(55) 



On the other hand, replacing k by ii' in H38|) (large K 
approximation) and M(0) by in Hij9|l . the mode con- 
dition (|22|l gives 



u_l(0) LOol^^nis Klsi + A 



(56) 



Equating ()54|) to H56() gives an estimate of the spin wave 
mode frequency u) = ojr -f We find 



UJl 



awR 



7'PA' 



Q ZZ 

30 



■A 



12 



Q ZZ 

jo 



1 



(57) 



The second line in H57|l confirms our expectation about 
the behavior of the mode frequencies and thus rational- 
izes the topology of Figure 0| Moreover, this analytic 
formula for lob. agrees quantitatively with the real part of 
the mode frequencies we find from our numerical solution 
of the full equations of motion. Unfortunately, l|57l) also 
predicts ui > 0. That is, (|53|l predicts only stable spin 
waves. These modes are unstable in the full solution of 
(|52|l . but only barely so (uji is negative but very small 
in magnitude). In fact, the stability or instability of a 
given mode turns out to be an extremely sensitive func- 
tion of the phase of m^(0)/m_l(0) and the factor M'{z) in 
(|52|l . The latter is determined by the ferromagnetic spin 
flip length Zgf- We conclude that variation of the magne- 
tization with the longitudinal variable z is an essential 
ingredient for instability to occur in the single-interface 
geometry. 

Figure [S] shows the calculated field dependence of the 
critical frequency Wcrit, the critical transverse wave vector 
Kcrit, and the critical current density Jcrit, where preces- 
sion first occurs in the single interface case. Consistent 
with (|57|) . all three are finite as — > 0. Otherwise, Wcrit 
and Jcrit vary roughly linearly with field, as they do in 
experiments performed on single films, ^ until the instabil- 
ity shuts off for large enough H. The order of magnitude 
of the critical current (« 10^ A/cm^) is also consistent 
with the single film experiments (and roughly ten times 
larger than the critical current seen in the trilayer pil- 
lars studied in Ref. Note, finally, that LUci-it deviates 
systematically from the Larmor frequency jH. 
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FIG. 6: Current dependence of the spatial variation of the 
magnetization for the three lowest spin wave modes of the 
ferromagnetic thin film shown in the middle panel of Figure^ 
Sohd curves are Mx{z). Dashed curves are My{z). There is 
almost no mode mixing for Jo = 0.05 X 10® A/cm^ 




ot:^ . . ^ . . . ^ . I 

0.00 0.05 0.10 

H(T) 

FIG. 5: Dependence on external magnetic field H of the criti- 
cal precession frequency Uait divided by the Larmor frequency 
(top panel); the critical parallel wave vector -ft'crit (middle 
panel); and the critical current density Jcrit (bottom panel) 
for the single interface geometry of Figure ^a). 



B. Thin Film Geometry 

We remarked at the end of Section IIIII that the in- 
terfacial spin torque associated with current flow from a 
ferromagnet to a non-magnet is opposite in direction to 
the interfacial spin torque associated with current flow 
from a non-magnet to a ferromagnet. As pointed out by 
Polianski and Brouwerji^ this implies that zero net spin 
torque acts on a uniformly magnetized film attached to 
identical leads at identical interfaces. Those authors used 
leads of different lengths to break the mirror symmetry 
of this geometry. We discuss that case in more detail in 
the next section. 

In our calculations for a ferromagnetic film sandwiched 
between two semi-infinite, non-magnetic leads (middle 
panel of Figure^, the natural spatial variation of the 
magnetization of the spin wave modes breaks the longi- 
tudinal symmetry. Additional spatial variation arises be- 
cause there is current-induced mixing amongst the modes 
(FigureEJ. As in the single- interface problem, the spatial 
variation M{z) is essential for precession if we assume a 



symmetrical geometry where the scattering properties of 
the two interfaces are identical. We do this here, even 
though the formulae given in Section IVll allow for a more 
general, asymmetrical situation. 

Figure [7| illustrates precessional instability diagrams 
for a ferromagnetic film of thickness t — AO nm. For 
each value of external magnetic field, the integer that 
labels each tear drop shaped region indicates the mode 
that goes unstable in that region. Notice that the area 
of the envelope in the Jq — K plane that encloses all the 
H = 0.05 T teardrops is much larger than the area of 
the H = 0.05 T instability region for the single-interface 
problem sketched in Figure^] The spin transfer torque is 
roughly the same in the two cases, but it is much easier 
to excite spin waves in the film because there is much 
less ferromagnetic "mass" present to contribute to the 
damping. 

The field dependence of the critical current is 
shown in Figure |H1 for several film thicknesses. A rela- 
tively large current density is needed to excite the 10 nm 
film because the symmetry breaking required to gener- 
ate a net torque occurs only when the lowest frequency 
spin wave mode (which has uniform magnetization in the 
z-direction) is perturbed by current-induced mixing with 
the next-highest spin wave mode (which has non-uniform 
magnetization in the z-direction). This is shown by the 
leftmost panel of Figure IHl We will use H55I) as the char- 
acteristic frequency associated with the driving current 
and (|51|1 to estimate the frequency splitting between the 
n = 1 and n — 2 film modes. Setting these equal, we 
get a field-independent estimate of Jcrit of the order of 
10^ - lOi" A/in. 

The critical current density decreases in Figure [S] as 
the film thickness increases to 20 nm and then to 40 nm 
because the mode splitting decreases in typical particle- 
in-a-box fashion. We will refer to this below as a quantum 
size effect. Then, beginning with the 40 nm film, a new 
phenomenon occurs. For large enough H , the perturbed 
n = 1 precession mode never goes unstable (there is no 
tear drop labelled "1" in the 77 = 2 T panel of Figure[7J). 
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FIG. 7: Instability diagram for a t = 40 nm ferromagnetic 
film in the geometry of Figure^b) for three values of external 
magnetic field. The integers label the thin film precessional 
mode that goes unstable in each region. 



FIG. 8: Critical current as a function of field for several 
thicknesses of the ferromagnetic film in the geometry of Fig- 
ure 0^b). The non-smooth behavior of the curves for the 
thicker films is associated with transitions between instabili- 
ties associated with different thin film modes. 



Instead, an upward tick in Jcrit occurs when the per- 
turbed n — 2 spin wave mode begins to precess. For the 
80 nm film, this phenomenon moves to lower values of 
H and successively higher spin wave modes become the 
first to become unstable. This transition of instability 
from one mode to the next highest is responsible for the 
jagged appearance of the Jcrit{H) curves for the thicker 
films. 

The nearly vertical line segment in Figure |S1 is the 
t = oo single-interface result redrawn from the lowest 
panel of Figure It is smooth because the region of 
instability for each mode becomes increasing small (as a 
function of H) until, in the limit of a semi-infinite film, 
every value of H leads to the linear instability of a differ- 
ent precession mode of the film. This curve terminates 
at a critical field i?crit(^ = oo). The 120 nm film curve 
shows that i?crit(0 increases as the film thickness de- 
creases from t — oo. 

Finally, FigureElsummarizes the thickness dependence 
of Jcrit for three values of external magnetic field. For 
H = 0.1 T, the critical current decreases monotonically 
as t increases as the quantum size effect argument sug- 
gests. However, to reach the t oo limit (horizontal 
arrow), Jcrit must eventually increase. As the iJ = 1 T 
and H = 3T curves show, this occurs in jumps each time 
a mode becomes absolutely stable and the next highest 
mode initiates the precession. Jcrit (^) oscillates in this 
regime because the quantum size effect operates as long 
as one particular mode is responsible for the instability. 

The iJ = IT and JT = 3T curves in Figure H termi- 
nate at a critical thickness tcrit{H) that decreases with 
increasing magnetic field. This is consistent with the fact 
that the semi-infinite ferromagnet only supports current- 
induced precession for relatively small values of external 
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FIG. 9: Critical current as a function of ferromagnetic film 
thickness for three values of external magnetic field in the 
geometry of FigureQ^b) . A jump occurs whenever a new mode 
is responsible for the precession. The horizontal arrow is the 
H = 0.1 T critical current for a semi-infinite ferromagnetic 
film. 

magnetic field. 

VIII. EFFECT OF FINITE LEADS 

Polianski and Brouwer«i^ studied current-induced mag- 
netization precession for the geometry depicted in the 
lowest panel of Figure ^ The unequal length leads 
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that break the symmetry of the problem are each at- 
tached to an electron "reservoir" . By definition, the elec- 
trochemical potential takes its equilibrium value inside 
each reservoir. Therefore, the spin accumulation (a non- 
equilibrium property) must go smoothly to zero when we 
enter each reservoir from the adjacent lead. These au- 
thors also fixed the magnetization of the ferromagnetic 
film to be uniform in the z direction. 

In this section, we study this reservoir geometry us- 
ing our formalism, both with and without the additional 
restriction to uniform longitudinal magnetization in the 
fcrromagnet. Our expectation is that the approximations 
of Ref.[l2 will be adequate if the ferromagnetic film thick- 
ness t is less than the exchange length ?ox- 

Let the non-magnetic leads have lengths t'^ and as 
shown in Figure ^c). It is straightforward to check that 
the solution of (|11() in, say, the right lead that satisfies 
the reservoir boundary condition of zero longitudinal and 
transverse spin accumulation at z = i -I- gives at the 
interface 

m,{t) = m+ (cxp[-<VZ,f] - exp[tV/,f]) 

mj.(t) = mi (exp[-Kt^] - exp[Kt^]) . (58) 
The corresponding spin current is 

Q..{t) = (i?m+/;,f)(exphtVZ,f]-f exp[t^//sf]) 

Q^M = (kI?to|) (exp[-Kt^] +exp[Kt^]) . (59) 

These can be recast as a relation between the spin accu- 
mulation and the spin current at the interface: 

QzyXt) = -^m,{t) COth{t^/kf) 
Q^^it) = -KDm^{t)coth{Kt^). 



(60) 



These formulae are the same as those used previously ex- 
cept for multiplicative factors depending on the length of 
the leads. We recover (|32|l and H3t)|) in the limit ^ oo. 
Very similar results apply to the left lead and interface. 

The modified spin currents (|60|) can be inserted into 
the solution of the longitudinal spin accumulation and 
then the boundary conditions for the interfacial spin 
torque. Wherever D appears explicitly (i.e. Igf and n 
are not affected) it can be replaced by D times the factor 
given above. Otherwise, the full calculation for the reser- 
voir geometry proceeds as outlined in previous sections. 
Figure [TUI illustrates results for two films with = oo. 
The = 450 nm is weakly lead asymmetric in the sense 
that the factor coth(t^/^sf) in H60() takes the value 1.3. 
The = 45 nm film is strongly lead asymmetric because 
coth(i^/^sf) has the value 10. 

From relatively small values near t — (see below), 
Jcrit in Figure ^1 increases initially as the film thick- 
ens due to increased damping from the additional ferro- 
magnetic material. This trend continues until the low- 
est mode of uniform magnetization becomes absolutely 
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FIG. 10: Critical current as a function of ferromagnetic film 
thickness in the asymmetrical geometry of Figure H^c) with 
= oo. The dashed curves are a uniform magnetization 
approximation. The solid curves are the resufis of the full 
calculation. The breaks in the solid curves correspond to the 
onset of precession from successively higher film modes. The 
external field is H = 3 T. 



stable and the next higher spin wave mode initiates pre- 
cession. Thereafter, the t dependence of the critical cur- 
rent density is determined by a competition between the 
quantum size effect (which tends to decrease Jcrit at t in- 
creases) and increased damping (which tends to increase 
Jcrit as t increases). 

A comparison of Figure El with Figure ^| shows that 
the very thin film limit of Jcrit depends very sensitively 
on the asymmetry of the leads. When the leads have 
equal length (FigureOJ, the spin current and torque come 
exclusively from the DzM'{z) term in H41|l . This implies 
a large critical current density because a large amount of 
energy is needed to mix the modes when i — > [see (l5T|l ]. 
By contrast, the critical current density is relatively small 
when the leads have unequal length (Figure [Tn|) because 
(even if the magnetization is uniform) there is incomplete 
cancellation between the (oppositely signed) spin transfer 
torques generated at the two ends of the ferromagnetii^ 
Figure lTTl shows that Jcrit actually increases exponentially 
with t^ when t^ = oo. 

We now consider an approximation to these full calcu- 
lations that forces the magnetization in the ferromagnetic 
film to be uniform in the z direction. This occurs auto- 
matically in our numerical work if the exchange length 
Zox and the spin diffusion length 4f are made much larger 
than any thickness t of interest. The results are shown as 
dashed curves in Figure^l As expected, the full and ap- 
proximate calculations agree quantitatively only when t 
is smaller than or approximately equal to the true Co ex- 
change length of 3.9 nm (see TablePl. The agreement for 
t = 10 nm is somewhat accidental because (as the break 
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FIG. 11: Log-linear plot of critical current in the geometry of 
Figure P^c) as a function of when — oo. The thickness 
of the ferromagnetic film is t = 2.5 nm. 

in the curve indicates) the precession at this thickness is 
not associated with a state of uniform magnetization in 
the longitudinal direction. 

On the other hand, FigurelTTilshows that for our choice 
of parameters, the uniform magnetization approximation 
gets relatively better (for t » lex) as the lengths of the 
two leads become more dissimilar. Qualitatively, this oc- 
curs because there are two sources of net spin transfer 
torque in this problem: the longitudinal asymmetry of 
the leads and the longitudinal asymmetry of the magne- 
tization. If the first dominates the second (strong lead 
asymmetry), it does not really matter which mode goes 
unstable and the uniform magnetization approximation 
will be acceptable. If the second dominates the second 
(weak lead asymmetry), the mode identity is crucial and 
the approximate calculation fails badly. 



IX. LIMITATIONS & FUTURE DIRECTIONS 

In experiment, the large current densities required to 
produce current-induced precession are achieved using 
point-contact, nanowire, or nanopillar geometries where 
conductor dimensions in the direction perpendicular to 
the current flow are "necked-down" from micron-scale 
dimensions to hundred nanometer-scale dimensions. The 
reservoir geometry of Figure ^c) is an approximate way 
to model this behavior within the context of a strictly 
one-dimensional calculation. Nonetheless, it remains true 
that wc have ignored the fact that the current only acts 
over a small area. The small area is either a small part 
of a large film, or a small finite element. In both cases, 
we have ignored the behavior of the current at the edges. 

For a point contact geometry, we have ignored how 
the current spreads into the film, how the variation in 



the magnetization decays into the film both laterally and 
in the interface direction, and how these interact. For 
a nanowire or nanopillar geometry, we have ignored the 
behavior of the current at the edges. Clearly, these three- 
dimensional effects will need to be considered if calcu- 
lations of the sort presented here are ever to be truly 
predictive in the sense of device design. 

Our calculations take no account of magnetostatics. 
The necessary computations are a bit complex because 
the magnetic self-energy is a non-local function of the 
spin density. It is true that the dominant effect can be 
modelled as magnetic anisotropy (which we have also ig- 
nored), but the known effects of magnetostatics on spin- 
wave dispersion suggest the possibility of other effects for 
the present problem. 

More significantly, we have focused exclusively on 
the initial, linear instability of a current-carrying nano- 
magnet toward precession. In reality, non-linear effects 
cut off the exponential growth of the mode amplitudes 
we find. This means that the final magnetic dynamical 
state may have little in common with simple precession, 
particularly for magnets whose size significantly exceeds 
an exchange length. To our knowledge, little or nothing 
is known about spatially inhomogeneous magnetic states 
of this type in the present context. 

Finally, we have performed all of our calculations at 
zero temperature. An extension to finite temperature 
clearly would be interesting in light of recent experiments 
on thermally activated magnetic reversal induced by a 
spin-polarized current^ 



X. SUMMARY 

We have combined a drift-diffusion description of 
transport with a Landau-Lifshitz-Gilbert description of 
magnetization dynamics to demonstrate that unpolar- 
ized current flow from a non-magnet into a ferromagnet 
can produce a precession-type instability of the magne- 
tization. Neglecting magnetic anisotropy and magneto- 
statics, our numerical calculations show (for both single- 
interface and thin film geometries) that a linear insta- 
bility occurs when both the current density and the ex- 
citation wave vector parallel to the interface are neither 
too small nor too large. The critical current density gen- 
erally increases with increasing magnetic field until the 
instability is quenched. Compared to the single-interface 
geometry, the film geometry supports a precessional in- 
stability over a much larger range of external magnetic 
field. In general, a particular spin wave mode initiates 
the instability when the mode just below it (in energy) 
becomes absolutely stable against current-induced pre- 
cession. 

In a lead-film-lead geometry, it is necessary to break 
mirror symmetry to generate precession because equal 
and opposite spin transfer torques are generated at the 
two ends of the ferromagnetic film. In our model, the 
required symmetry breaking is provided by spin wave 
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modes of the film whose magnetization varies along the 
direction of current flow. If the magnetization is artifi- 
cially constrained to be longitudinally uniform, asymme- 
try must be introduced in the length or properties of the 
leads or in the scattering properties of the two interfaces. 
When the spin wave mode asymmetry generates the spin- 
transfer torque, the critical current density first decreases 
(a quantum size effect) and then slowly increases (as suc- 
cessive modes become stable against precession) as the 
film thickness increases. When interface or lead asymme- 
try dominates the spin transfer torque, the critical cur- 
rent is small (compared to the mode-dominated regime) 
and slowly increases due to increased damping as the film 
thickness increases. 
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Appendix B: Transport Boundary Conditions 

The interfacial boundary conditions 123() for the drift- 
diffusion equations l(2Jl, (QJ, 0, and (O can be derived 
most convincingly by averaging the boundary conditions 
for the semi-classical distribution functions (/|(k, r) and 
(7l(k, r) used in a Boltzmann equation treatment of spin 
transport .-^^ A simpler alternative is to proceed phe- 
nomcnologically and write kinetic equations that connect 
the currents to the densities. Temporarily, we assume 
that the ferromagnetic magnetization M = Mu is paral- 
lel to the interface normal, i.e., parallel to the z-axis. 

In linear approximation, the longitudinal particle cur- 
rent that flows through the interface for each spin type 
is 



(67) 



31 - h{^,'^^-^,J^) = lk,{a,5n^^-a,5nl^) 



Appendix A: Exchange Magnetic Field 



The total exchange energy of a semi-infinite ferromag- 
net that extends from z = zl to z = zr, is 

Sex - ^ / <PR / dz{do,M.){doM). (61) 

The quantity rus occurs here because M denotes spin 
density rather than magnetization. Focus on the a = 
z component and let M — > M + 77. This induces the 
variation 



2A, 



d^R dz{d,_'n)-{dM)- (62) 



Integrating by parts gives 



SB,,, = 



-2A, 



d'R 



-2 A, 



dzri-dlM.-['n-dM)] 



■J d^R dzr]-{dlM. 



(63) 



+ dM[5{z-zi^) - 5{z-z^)]}. (64) 

The exchange field iJcx = -7"^((5i;cx/'5M). Therefore, 
since J7 = 5M, 

2^ox 



7m| 



{dim + dM W - ^l) - 6{z - ZR,)]} . 



(65) 

The a = X and a = y components of H61|l do not generate 
surface terms so (|65|) generalizes to 

2A„ 



7m~ 

2j4ex 

7to2 



[6iz - zl) - S{z - zr)] OM- (66) 



Here, fc| and fc| are kinetic coefficients, n^^, 



,,NM ,,FM 



and /i™ are electrochemical potentials, and oo, a|, and 
ai are twice the inverse density of states per spin at the 
Fermi level for electrons in the non-magnet and for ma- 
jority and minority electrons in the ferromagnet, respec- 
tively. Using H67() , the same algebra that led from JBJ to 
(0 now yields 

3 = 31 + 31 

= (fct + ki)aoSn^^ - (aifct + aiki)Sn™ 

2 2 

+ -^{k-i - ki)aom^ - -{a-ik^ - aiki)Sni. (68) 

The physical charge accumulation Suq at the interface 
can be partitioned arbitrarily between Sn^^ and Sri^^. 
If we choose 612 — Sn^^ and Sn^^ — 0, the boundary 
condition for the particle current becomes 



j = {ki^ + ki)ao6nQ + 2h ^{k^—ki)aor 
— 2h^^{a^k]- — aiki)Sm. 



(69) 



In exactly the same way, we derive the boundary condi- 
tion for the longitudinal spin current from 



JT -31 
i?i(fc| - 

{a^k^ 4 



ki)aoSn + (/s| + fc^)aoTO 
aiki)Sm. 



(70) 



Equations (|69|l and H7U|I reproduce the first and last 
boundary conditions in (|23() if we use the Boltzmann 
equation treatment of the boundary conditions^ to make 
the choice 



1-R. 



(27r)3?i 4R^ 



(71) 
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FIG. 12: The z-axis of unprimed coordinate system is per- 
pendicular to the material interface. The z-axis of the primed 
coordinate system is parallel to the magnetization M. 



We stated the boundary condition for the transverse 
spin current in Section III, leaving a constant of propor- 
tionality b unspecified in Hl()|) . As above, we choose this 
constant for consistency with the Boltzmann equation 
treatment. The result is 



QNM _ 



A 



FS 



aom_L, 



(72) 



where A-ps is the area of the Fermi surface projected onto 
the interface. The reflection coefficients in (|71|) do not ap- 
pear in (|72|) despite the fact that spins can rotate into the 
transverse direction upon quantum mechanical reflection 
or transmission at an interface. This is so because the 
rotation effect disappears after Fermi surface averaging,^ 
Tserkovnyak and co-workers2i discuss another contri- 
bution to due to "spin-pumping" from an adjacent 
ferromagnet. When the magnetization is time- harmonic 



with frequency w, the time derivative of the magneti- 
zation can have a component that is transverse to the 
magnetization and there is an additional term in H72|l 



Q 



NM 



A 



FS 



(27r)3?i 



aom_L 



Ays 

(27r)3 ~2 



U-L, 



(73) 



For our choice of parameters, this new term is numeri- 
cally much smaller than the other term except for high 
frequencies and low currents. Nevertheless, we have re- 
tained it for all the calculations presented here. 

As a final step, we must take explicit account of the fact 
that M is tilted slightly away from the interface normal 
z. This is not difficult if we define a a new coordinate 
system x'y'z' where z' points along M (Figure EJ. In 
that case, for example. 



NM 



A 



FS 



-aom_L' 



(74) 



(27r)3?i 

replaces (I72II as the boundary condition for the transverse 
spin current. However, because the direction of u_l in 
Figure IT21 varies along the interface when the excitation 
wave vector X 7^ 0, it behooves us to evaluate (|7^ using 
the original, interface-referenced coordinates xyz. 

The transformation law between the xyz and x'y' z' 
coordinate systems (Figure I12|) simplifies considerably 
when u\_ ^1, as is the case for our problem. Then, 



(75) 



for any vector A that lies close to the z axis. This rule 
shows that (|74|l is indeed identical to the transverse 
boundary conditions stated in H23|l . It also shows that 
the longitudinal boundary conditions stated above are 
unaffected by the fact that uj^ ^ 0. 





— Ax 


u±A^ 


Ay 


= Ay - 


uj_A:, 









1 E. B. Myers, D. C. Ralph, J. A. Katine, R. N. Louie, R. 

A. Buhrman, Science 285, 867 (1999); J. A. Katine, F. J. 

Albert, R. A. Buhrman, E. B. Myers, D. C. Ralph, Phys. 

Rev. Lett. 84, 3149 (2000). 
^ J. GroUier, V. Cros, A. Hamzic, J. M. George, H. Jaffres, 

A. Fert, G. Faini, J. Ben Youssef, H. Legall, Appl. Phys. 

Lett. 78, 3663 (2001). 
^ J.E. Wegrowe, D. Kelly, T. Truong, Ph. Guittienne, and 

J.-Ph. Ansermet, Europhys. Lett. 56, 748 (2001). 
" J. C. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996); 

J.C. Slonczewski, J. Magn. Magn. Mater. 195, L261 

(1999). 

^ L. Berger, Phys. Rev. B 54, 9353 (1996). 

® M. Tsoi, A. G. M. Jansen, J. Bass, W. G. Chiang, M. Seek, 

V. Tsoi, P. Wyder, Phys. Rev. Lett. 80, 4281 (1998); J. 

Z. Sun, J. Magn. Magn. Mater. 202, 157 (1999); J.E. We- 



grove, D. Kelly, P. Guitienne, Y. Jaccard, and J. P. Anser- 
met, Europhys. Lett. 45, 626 (1999). 

Y. Ji, G.L. Chien, and M.D. Stiles, Phys. Rev. Lett. 90, 
106601 (2003). 

* S.I. Kiselev, J.G. Sankey, I.N. Krivorotov, N.G. Em- 
ley, R.J. Schoelkopf, R.A. Buhrman, and D.C. Ralph 
cond-mat/0306259 W.H. Rippard, M.R. Pufall, S. Kaka, 
S.E. Russek, and T.J. Silva cond-mat/0309042 Previous 
experiments [M. Tsoi, A.G.M. Jansen, J. Bass, W.-G. 
Ghiang, V. Tsoi, and P. Wyder, Nature 406, 46 (2000)] 
demonstrated that external microwave irradiation induces 
new resistance anomalies in magnetic multilayers when the 
current exceeds a critical value. 

^ M.D. Stiles and A. Zangwill, Phys. Rev. 66, 014407 (2002). 

° M.D. Stiles and A. Zangwill, J. Appl. Phys. 91, 6812 
(2002). 



15 



Y. Tserkovnyak, A. Brataas, and G.E.W. Bauer, Phys. 
Rev. B 67, 140404 (2003). 

Ya. B. Bazaliy, B.A. Jones, and S.-C. Zhang, 
cond-mat/0009034 

J.Z. Sun, Phys. Rev. B 62, 570 (2000). 
^'^ Ya. B. Bazahy, B.A. Jones, and S.-C. Zhang, Phys. Rev. 
B 57, R3213 (1998). 

M.L. Pohanski and P.W. Brouwer, cond-mat/0304069 
A. Brataas, Y.V. Nazarov, and G.E.W. Bauer, Eur. Phys. 
J. B 22, 99 (2001). 

A paper that considers spin accumulation and spin current 
transverse to the ferromagnetic magnetization is A. Shpiro, 
P.M. Levy, and S. Zhang, Phys. Rev. B 67, 104430 (2003). 

1* A. Pert, J. Phys. C 2, 1784 (1969). 

X. Waintal, E.B. Myers, P.W. Brouwer, and D.C. Ralph, 
Phys. Rev. B 62, 12317 (2000). 

^° The first term on the right side of ^ shows that Q™ 
is positive deep in the bulk where Sm = 0. The diffusive 
correction can change its sign near the interface, but this 
does not occur for the physical range of parameters we use. 

'^^ See, for example, A.G. Gurevich and G.A. Melkov, Magne- 
tization Oscillations and Waves, (CRC Press, Boca Raton, 
1996), Chapter 7. 

J. Bass and W. P. Pratt, J. Magn. Magn. Mater. 200, 274 
(1999). 

Q. Yang, P. Holody, S.-F. Lee, L. L. Henry, R. Loloee, P. 



A. Schroeder, W. P. Pratt, Jr., and J. Bass , Phys. Rev. 
Lett. 72, 3274 (1994). 

L. Piraux S. Dubois, A. Pert, and L. Belliard, Eur. Phys. 
J. B 4, 413 (1998). 

N.W. Ashcroft and N.D. Mermin, Solid State Physics 
(Saunders, Philadelphia, 1976), Table 2.1. 

D. A. Papaconstantopoulos, Handbook of the Band Struc- 
tures of Elemental Solids (Plenum, New York, 1986). 

A range of values have been found experimentally, see for 
example: G. Shirane, V. J. Minkiewicz, and R. Nathans, J. 
App. Phys 39, 383 (1968) and Z. Frait, Phys. Stat. Solidi 
2, 1417 (1969). 

R. Pauthenet, J. Appl. Phys 53 8187 (1982). 

S. M. Bhagat and P. Lubitz, Phys Rev. B 10, 179 (1974). 

R.D.McMichael (private communication). 

E. B. Myers, F.J. Albert, J.C. Sankey, E. Bonet, R.A. 
Buhrman, and D.C. Ralph, Phys. Rev. Lett. 89, 196801 
(2002). 

M.D. Stiles (unpublished). 

Y. Tserkovnyak, A. Brataas, and G.E.W. Bauer, Phys. 
Rev. Lett. 88, 117601 (2002); A. Brataas, Y. Tserkovnyak, 
G.E.W. Bauer, and B.I. Halperin, Phys. Rev. B 66, 060404 
(2002); Y. Tserkovnyak, A. Brataas, and G.E.W Bauer, 
Phys. Rev. B 67, 140404(R) (2003). 



16 



